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Abstract 

Using the Mathematica program we calculate numerically the dif- 
ference of the diagonal matrix elements of the time dependent effective 
Hamiltonian for the neutral K meson complex. We consider the exactly 
solvable neutral K meson model based on the one-pole approximation 
for the mass density. The so-called Khalfin's Theorem is numerically 
examined. Some characteristic parameters for this system are also cal- 
culated. The results of all calculations are presented in the graphical 
form. The calculations are made assuming the total system is CPT - 
invariant and CP - noninvariant. 

PACS numbers: 03.65. Ca, ll.30.Er, 11. 10. St, 14.40.Aq 
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1 Introduction 

The neutral kaon system is probably one of the most interesting complexes 
of elementary particles. Using this system it was found in 1964 that the 
CP symmetry is violated pQ. The description of CP violation effects in this 
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system is based on the approximation proposed by Lee, Oehme and Yang 
(LOY) in [2]. This theory was then developed and intensively studied by 
Lee (see [3]) and by many other authors (see eg. [4]). Within the LOY 
approach, a non-hermitian Hamiltonian H\\ is used to study the properties 
of the particle-antiparticle unstable system [2, [HI SJ [5] 

#u = M-~r, (i) 

where 

M .\r . V I' (2) 

are (2 x 2) matrices acting in H\\, where H\\ is a two-dimensional subspace 
of the total Hilbert space of states H, spanned by the state vectors of K°, 
K° mesons. The M-matrix is called the mass matrix and T is the decay 
matrix. LOY derived their approximate effective Hamiltonian H\\ = H^qy 
by adapting the one-dimensional Weisskopf-Wigner (WW) method to the 
two-dimensional case corresponding to the neutral kaon system. Almost all 
properties of this system can be described by solving the Schrodinger-like 
equation 

ij^;*)|| = #ufoM)||, (t>to>-oo) (3) 

where we have used H = c = 1, and \if); t)\\ G TC\\. 

Within the LOY theory the physical states of neutral kaons are superpo- 
sitions of \K°) and \K°). They are the eigenvectors of Hloy, 

\K s )=p\K°) + q\K ), \K L )=p\K°)-q\K°) t (4) 

H\\\K S(L) ) = Hs{l)\K S {l))i (5) 

and correspond to short-living (the vector \K S )) and long-living (the vector 
\K L )) states of neutral kaons. We will use the following notations further on 
in this paper: \K°) = |1), \K°) = |2). 

One of standard results of the LOY approach is the following: In a CPT 
invariant system, i.e. when 

OH&- 1 = H, (6) 

(where = CPT, and H is the total self-adjoint Hamiltonian for the system 
containing neutral koans considered), there is 

= h L 22 ° y (?) 
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and 

MT = M 2 L 2 OY , (8) 

where: Mff Y = &(hM >Y ) and $l(z) denotes the real part of a complex 
number z ($S(z) is the imaginary part of z), and hJ-° Y = (S\Hloy\]) , (J = 
1,2). Another important prediction of the LOY theory is that the ratio 

k*) = K = 44^ = const and i^)i = i4i^i, (9) 

q 2 A 2 i{t) q 2 
when CP symmetry is violated, [CP,H] ^ 0, [3] - [8]. Here 

^(t) = (j|f/||(t)|k) ee (j|e-^l|k) (j, k = 1,2), (10) 

and U\\(t) is the evolution operator for the subspace 7i\\. 

The important result indicating some limitations of the LOY approach 
was obtained by Khalfin P El El Q33 QIl Q21- Khalfin found that in the 
exact theory there must be 

if r (t) = = const., then \r(t)\ = 1, (11) 



where 



A 3k (t) = (j\e- itH \k), (j, k = 1,2). (12) 



Result (fill ) is known in the literature as "Khalfin's Theorem". Using this 
result Khalfin hypothesized that beyond the LOY approximation one should 
expect new CP-violation effects [71 [8] and he tried to obtain some model es- 
timations of the possible magnitude of these effects. He found that the order 
of these effects should be 1CT 3 (see |8J). He obtained his estimation using the 
spectral language for the description of Ks,Kl and K°, K°, by introducing 
a hermitian Hamiltonian, H, with a continuous spectrum of decay products 
labeled by a, (3, etc., 

H\<f) a (m)) = m \(f) a (m)), (<j>p(m')\<f> a (m)) = S af3 5(m' - m). (13) 

Here H is the above mentioned total Hamiltonian for the system. H includes 
all interactions and has absolutely continuous spectrum. We have 

\K S ) = / dm V% !a (m)|i(m)), (14) 

JSpcc (ff) 
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and 



K L ) = dm y^cuz,, a (m)l^(m)) ; (15) 

JSpec (H) Q 

|j) = / dm y2u jia (m)\(/) a (m)), (16) 

./Spec (m 



where j = 1,2. Thus, the exact Ajk(t) can be written as the Fourier transform 
of the density pjk(m), (j, k — 1, 2), 



A ifc (t)= / dme-* m Vi*M, (17) 

J — oo 

where 

pjfc(m) = 2^a;* a (m)a; fcja (m). (18) 

The minimal mathematical requirement for Pjk{m) is the following: 
f_°° dm \pjk(m)\ < oo. Other requirements for Pjkijn) are determined by 
basic physical properties of the system. The main property is that the energy 
(i.e. the spectrum of if) should be bounded from below, Spec(H) = [m g , oo) 
and m g > — oo. 

Starting from densities Pjkijn) one can calculate Ajkit). In order to find 
these densities from relation ([TBI one should know the expansion coefficients 
u)j >a (m). Using physical states \Ks), \Kl) and relations ([4]) they can be ex- 
pressed in terms of the expansion coefficients c^s, a (m), ujs >a (m). Thus, as- 
suming the form of coefficients ujs )a {rn),u}s,a( m ) defining physical states of 
neutral kaons one can compute all Ajk(t), (j, k — 1, 2). 

The model considered by Khalfin is based on the assumption that (see 
formula (35) in |8j), 

sA m S ~ «f) I m ~ m S + l 2 

u L Am) = \——- ^— — p^, (20) 

V 111 \KsA m L ~ *-^)| m-m L + i^- 

where as$ and are the decay (transition) amplitudes and ^s(L)A m ) 
are, in general, some nonsingular "preparation functions". Khalfin found his 
above mentioned estimation choosing, for simplicity, the trivial form of the 
"preparation functions", £s(L),p( m ) = 1- 
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The discussion about the validity of the Khalfin's estimation of his new 
CP violation effect can be found in the literature (see eg. [12 CE2])- Our at- 
tention will be concentrated on the attempt to verify the size of the Khalfin's 
estimation performed in [TJ]. The calculation performed in [12] uses Khalfin's 
assumption that £s(L) )( g(m) = 1, strictly speaking, they use the assumption 
that in (USD, (EDI) there is 



T = 9{m - m g ) = [g(m - m g ) 



def f 1 if m > m g , 



Vg , 



(21) 



def f 

[0 if m < m 
Within this assumption one obtains, for example, that 

. r+oo 

A ss (t) = (K s \e- ltH \K s ) = / dm p ss (m) e' ltm , (22) 

J — CO 

where 

Tg S 

Pss(m) = g(m-m g ) ^ — , (23) 

(m — mg) 2 + 

S = J2MKs^a)\ 2 , (24) 

a 

and so on. 

The form of density pgg(m) defined by (l23l) is not the most general one. 
In more realistic models functions Uj ta ( m ) an d of type LOs,p{ m )^ u;L ,p( m ) ^ ea ^ 
to the densities p{m) which general form is similar to (|23l) with g{m — m g ) 
and S = 5(m) having more involved form [T3l [i~4] . In the general case the 
threshold factor g(m — m g ) describes the behavior of p(m) for small m, (i.e., 
for m ~ m g ) and it is responsible for the long time properties of amplitudes 
of type Ajk(t) and Ass{t)- The second factor in the formulae of type (1231) 
having the Breit-Wigner form results from the pole structure of functions 
of type ujs,i3(ni),ujL 1 i3('m) (see ( 1T91) . (1201) ) defining densities p(m) and it is 
responsible for the form of Ajk(t),Ass{t) etc. for the intermediate times 
(i.e. it is responsible for the exponential part of the survival probabilities). 
The third factor, i.e. the factor corresponding to S(m) ensures the suitable 
behavior of p{m) for m — > oo. 

For simplicity, it is assumed in [12] that m g = 0. So all integrals of type 
( 1221) and ( 1171) are taken between the limits m = and m = +oo. All these 
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assumptions made it possible to express amplitudes of type Aj k (t) in p2] 
in terms of known special functions. The same assumptions were used in 
[T5] (see [H], relations (37) - (39) and (42) - (47)) and will be used in this 
paper. Note that putting g{m — m g ) = 1 in (1221) leads to strictly exponential 
form of amplitudes of type Assif) as functions of time t. On the other 
hand, keeping g(m) in the assumed simplest physically admissible form ( 1211 
results in the presence of additional nonoscillatory terms in amplitudes of 
type Ass(t), ALbif) etc. and thus in amplitudes Ajk(t) as well (see [T2l fT5] ). 

In [15] the analytical formulae for Ajkit) obtained in [12] were used as 
the starting point to find analytical expressions for matrix elements of the 
effective Hamiltonian for this model for t — tl and then to obtain a numerical 
value for the possible consequence of the Khalfin's Theorem analyzed in [16J. 
It is found there that, contrary to the standard LOY result (JJJ), the diagonal 
matrix elements of the exact effective Hamiltonian for neutral meson complex 
cannot be equal if CPT symmetry holds but CP symmetry is violated. We 
found in [15] that 

K(/in(t ~ r L ) - h 22 (t ~ t l )) ~ -4.771 x 10~ 18 MeV, (25) 

9?(/iii(* ~ t l ) - h 22 {t ~ t l )) ~ 7.283 x lO'^MeV (26) 



and 



|3?(/in(t ~ t l ) - h 22 (t ~ t l ))\ _ m K o - m^ 2J 



TFlaverage ^Haverage 



10" 21 , (27) 



where hjk(t) = (j\H\\(t)\k) ij,k = 1,2) and H\\(t) is the effective Hamilto- 
nian. These results were obtained analytically for the considered model for 
the neutral kaon system in the case when the total system is CPT - invariant 
but CP - non-invariant (equations (68), (69) and (70) in [IH]). The estima- 
tions (|25|) - (l27l) were obtained by inserting (|20l) and related rris — nriL — 
497.648MeF, Am = 3.489 x lQ- 12 MeV, t s = 0.8935 x 10^ 10 s, 



m, 



average 



TL = 5.17 x 10~ 8 s, 7 L = 1.3 x 10- u MeV, 75 = 7.4 x 10~ 12 MeV in formulae 
of type ([T9l . In this paper we will use the same experimental data. We will 
also use the same notations and definitions as in |15|: 



T s T L 

7s = 1l = -y, Am = m L - m s , (28) 



and so on. Note that results (12511 - (1271) agree with the general result obtained 
in \W\. 
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The detailed analysis of the matrix elements of the effective Hamilto- 
nian for the K° — K° system shows that the non - zero difference between 
the diagonal matrix elements of the effective Hamiltonian in the consid- 
ered model is caused by the nonzero contribution into pjk(m), (fT8l) . com- 
ing from expressions for (Ks\e~^^ \Kl) and (Ki\e~ * \Kg) and by the 
non-oscillatory terms in the formulae for the amplitudes of type (fT7l) for 

transitions: K° < ► K°, K° < ► K° K° < ► K°. It is not difficult to verify 

that neglecting the mentioned nonzero contribution and dropping all these 
non-oscillatory terms leads to the zero difference of the diagonal matrix ele- 
ments of the effective Hamiltonian in the considered case. This is because, in 
fact, dropping these non-oscillatory terms is equivalent to replacing in f fT7l ) 
densities Pjk{m) defined by (fT8l ) - ( 1211 with densities defined by the new 
function gww( m ) instead of g(m) given by (l2Tj) such that gww{ m ) — 1 f° r 
all — oo < m < +00. Thus, the integrals, e.g. in formulae of type (l22l . are 
taken between the limits m = —00 and m = +00 with densities of type ff23l) 
having the Breit-Wigner form (and not truncated for m < m g = 0) which 
leads to strictly exponential form of, e.g. |As,s(t)| 2 and the like. The effective 
Hamiltonian H\\ (£) obtained in such a case is the LOY effective Hamiltonian, 
Hloy- 

In this paper we continue searching for the properties of the model ana- 
lyzed in [15]. The aim is to show how the difference of the diagonal matrix 
elements (hn(t) — ^22^)) discussed in p[5] and some other parameters de- 
scribing neutral kaons (including r(t), f fTTl) ) change in time t in the case of 
preserved CPT and violated CP symmetries. The paper is organized as fol- 
lows. In Section 2 we collect the formulae for the matrix elements of the 
effective Hamiltonian necessary for further analysis. Section 3 contains a 
numerical verification of the Khalfin's Theorem (fTTTl . We show there how 
this Theorem "acts". In Section 4 the value of the difference of the diagonal 
elements (h n (i) — h 2 2(t)) is shown as calculated at t = r L with the use of 
the Mathematica. Also, in this Section the time dependence of the real and 
imaginary parts of the diagonal matrix elements of the effective Hamiltonian 
(hn(t) — h 2 2(t)) in graphical form is given. In Section 5 the eigenvalues of 
the effective Hamiltonian fiL^t), /is(t) are calculated and the parameters of 
the violation of the CP symmetry EL(t),es{t) are estimated. We also show 
graphically the time dependence of all the calculated quantities there. In 
Section 6 we check the correctness of our results by verifying the relation 
(^i(t) + /is(£) = hn(t) + h 2 2(t)) known from the literature. Section 7 con- 
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tains a discussion of the results obtained in Sections 3-6 and some remarks 
concerning the experiments with neutral kaons. 



2 Matrix elements of the effective Hamiltonian 

General conclusions concerning properties of matrix elements of the effective 
Hamiltonian H\\ (£) can be drawn using the following identity [H] 



HM)=i^j&[A{t))-\ (29) 

where all matrix elements Aj k (t), (j,k = 1,2) of the matrix A(t) can be 
calculated, e.g. by means of ( |T7l ). 

Using ff29l . one can calculate all the matrix elements of the effective 
Hamiltonian H\\ which may now be written as 

M«> = ^TTV ' ' " ^ ' *.<*)), (30) 



MO = T777T77T ' ( ' -M<) + ■ A u (t) ) , (31) 



i 


(dA n (t) 


detA{t) 


v at 


1 .| 


r ^n(t) 


detA(t) 




z 


f&4 21 (£) 


detA(t) 


V 0* 




r 9^21 (*) 


detA(t) 


v at 



where 



M*) = "jTTTTIV ' ( • ^(t) - • A 21 (t) ) , (32) 



M*) = TT^TTTT ■ ( ■ + ■ ^ii(t) ) , (33) 



detA(t) = A n (t) ■ A 22 (t) - A 12 (t) ■ A 21 (t) . (34) 



So, having these relations and inserting analytical expressions for Ajk(t), 
ffTTl) . calculated in [15j within the assumptions ([6]) and [CP,H] ^ 0, one 
obtains all matrix elements hjk(t) for the model considered. Next, such 
obtained analytical formulae for hjf.{t) can be used for numerical calculations 
of some parameters characterizing neutral kaons for instants of time changing 
in given time intervals. 
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3 Numerical examination 
of the Khalfin's Theorem 



It seems to be interesting to verify how the Khalfin's Theorem ( fTTl) acts in 
the system on neutral mesons. To see this, we can use amplitudes A 12 (t) 
and A 2 i(t) calculated within the model considered in pj)] in the case of con- 
served CPT and violated CP symmetries. It is not difficult to calculate the 
modulus of the ratio using numerical methods. The results of such cal- 
culations are presented below in FigH] and FigfJl There are y(x) = 
and x = ■ t in these figures. 



1.0 



0.9 




Figure 1: The time dependence of the absolute value of y(x) = \r(t)\ 
in x G (0.01, 103). Here and in all other Figures: x = ^ • t. 



\ A 12 (t) 
\A 2 i(t) 
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1.0 



0.9 



4 ri l ^ NfH^irt 



0.Z 0.4 0.5 0.8 



Figure 2: The time dependence of the absolute value of y(x) = \r{t)\ = | 
in x e (0.1, 1). 



These figures show that if one is able to measure the modulus of the ratio 
only up to the accuracy of the order 10 -15 then one sees this ratio as 
a constant function of time: for x G (0.01, 103) we find that 

ymax{x) - Vmin{x) = 3.33067 x lO" 16 , (35) 

where 



| max > 



Vmaxi^x) 1^(^)1? 
Vmin\ x ) = \ r if)\min- (36) 



4 The difference of the diagonal matrix ele- 
ments {hn{t) - h 2 2{t)) 

Assuming that the CPT symmetry is conserved in the system under consid- 
erations ([CPT,H] = 0) and using the necessary relations from [T5l fT2l [T6] 
one finds the general form of the difference of the diagonal matrix elements 
of the effective Hamiltonian. It has the following form 
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h n (t) - h 22 (t) 



X(t) 



(37) 



detA(t) 



where 



X(t) 



= i ■ 



dA 21 (t) 
dt 



■A 12 (t) 



dA 12 (t) 
dt 



A 21 (t) 



(38) 



and detA(t) is defined in (|34l) . 

Our analytical result in the one-pole approximation obtained in [15] for 
t = tl can be written as 



Mtl) - h 22 (r L ) ~ (-4.771 x 1(T 18 + z • 7.283 x l(T 16 )MeV. (39) 



The numerical result for t = tl in the one-pole approximation obtained 
using the Mathematica has the following form 



hu{r L ) ~ h 22 (T L ) ~ (-7.129 x 10" 17 + % ■ 1.986 x lO~ n )MeV. (40) 



It is seen, that the difference between (|39l and (1401) is small and it may be 
attributed to finite accuracy of numerical calculations performed by Mathe- 
matica. No approximations have been used in the analytical calculations. 

Putting Ajkit) (j,k = 1,2) given by f fT7l) into ( 1371) and using the energy 
density Pjk{m) found in [15], the difference (hn(t) — h 22 (t)) can be calculated 
as the function ot time t. The results of our calculations are presented in a 
graphical form. The figures below show the time dependence of the real and 
imaginary parts of the diagonal matrix elements of the effective Hamiltonian 



(hn(t)-h 22 (t)). 
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Figure 3: The time dependence of the real part of the diagonal matrix el- 
ements of the effective Hamiltonian y = Re(hu(x) — h 2 2(xj) in the range 
x G (0.001, 1). 



Note that in FigH we have Re(h n (t) - h 22 (t)) = Re(h n (r L ) - h 22 (r L )) 
for x = 1. 



12 



Y 




x 



-2 • 10 



-19 



-4 -10 



-19 



-6-10 



-19 



-8 • 10 



-19 



-1 -10 



-18 



Figure 4: The time dependence of the real part of the diagonal matrix ele- 
ments of the effective Hamiltonian y = Re(h n (x) — h 2 2(x)) in the interval 

x e (1,103). 

At this point it should be explained that a more accurate analysis of the 
results of the calculations which lead to Figj3] and FigfJ] and the use of a 
larger scale show that the obtained curves are not so smooth as can be seen 
in Figj5] and FigJB] but they are similar to curves in Figj5] and FigJBl 
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Figure 5: The time dependence of the imaginary part of the diagonal matrix 
elements of the effective Hamiltonian y = Im(hn(x) — h 2 2(x)) in the interval 
x e (0.001,1). 
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Figure 6: The time dependence of the imaginary part of the diagonal matrix 
elements of the effective Hamiltonian y = Im(hn(x) — h 2 2(x)) in the interval 
x E (1,103). 
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5 Calculation of ^(t), fJ>s(t) an d £ L(t),£s{t) 

The eigenvalues of the effective Hamiltonian //l(£), fis(t) can be written as 

m 

tiL{t) = Kit) - h{t), (41) 

fi S (t) = h (t) + h(t), (42) 

where 

/i (t) = i-(Mt) + M*)), ( 43 ) 

h(t) = y/hl(t) + h 12 (t)-h 21 (t) (44) 

and 

h z {t) = y(hu(t)-h 22 {t)). (45) 
From (IHD and ((42j) we have 

M*) + fJ. L {t) = h n (t) + h 22 (t) = Tr(H\\(t)). (46) 

Relation (1461) does not depend on any approximations and it is always true 
for every (2 x 2) matrix. Inserting (l30l) - (l33l) into (14B and (1341) and then 
using ( fTTl) and performing all integrations of type (fl~7l) one can obtain for 

t = T L 

Hl{j l ) ^ (497.648 - i ■ 4.458 x l(T 13 )MeV, (47) 

and 

/i 5 (r L ) ~ (497.648 - i ■ 2.471 x l(T 13 )MeV. (48) 

The general formula for fj,L(s)(t) can also be written as follows 

i 

HL(s){t) = m L{S ){t) - ^ 7l(S)(*)- (49) 
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The results of our calculations for the real part and imaginary part in 
( 1471) and (1481 ) are rounded to the third decimal place. It should be noted 
that the real part in (1471) and the real part in ([48]) differ in the fourteenth 
decimal place. The above mentioned result corresponds with the fact, that 
m s 7^ f^L and there is \rriL — ms\ ~ |7s|, [201 E]- 

The time dependence of and /^l(^) is given below. 
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Figure 7: The time dependence of the real part of : y = i?e(^i(a;)) in 

the interval x G (0.001,103). 
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Figure 8: The time dependence of the real part of Hs{ x ) '■ V = Re{^s{ x )) m 
the interval x e (0.001, 103). 



Expansion of scale in Figd and FigJS] shows that continuous fluctuations 
with amplitudes of the order of 10~ 14 appear. 
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Figure 9: The time dependence of the imaginary part of /j-l{x) 
Im{nL{x)) in the interval x G (0.001, 103). 
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Figure 10: The time dependence of the imaginary part of fJ>s(x) : y = 
Im(fj,s(x)) in the interval x e (0.001, 103). 
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We have the following formulae (see, e.g. [18J) 

h 2 i(t) - h 22 (t) +fi L (t) 



h 2 i(t) + h 22 (t) -fi L (t) 
h 2 i(t) +h 22 (t) -ns(t) 



h 2 i(t) -h 22 (t) +ns(t) 
and we get for t = r L 

e L (r L ) ~ -1.0000000000184743' + i ■ 0.0' 

and 

s s (t l ) ~ 1.0000157759810688' - i ■ 0.0' 
The time dependence of Eiif) and esif) are presented below. 

y 
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Figure 11: The time dependence of the real part of el(x) : y = i?e(e^( 
the interval x G (0.001,103). 
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Figure 12: The time dependence of the real part of es(x) : y = Re(es(x)) in 
the interval x E (0.001, 103). 

Expansion of scale on FigfTI] and Fig{T2] shows, that continuous fluctua- 
tions with amplitudes of the order of 10~ 12 appear here. 
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Figure 13: The time dependence of the imaginary part of Sl{x) : y 
Im(e L (x)) in the interval x £ (0.001, 103) 
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Figure 14: The time dependence of the imaginary part of es(x) : y = 
Im(e s (x)) in the interval x G (0.001, 103) 
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From the formula 

e(t) = ±>(e L (t)+e s (t)) (54) 

we have for t — tl 

e(r L ) ~ 7.888 x 1(T 6 - i • 0.0'. (55) 
The absolute value of s(tl) 

\e(t l )\ ~ 7.888 x 10~ 6 . (56) 
The figures below present the time dependence of the absolute value of 
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Figure 15: The time dependence of the absolute value of e(t): y = \s(t)\ in 
the interval x e (0.001, 103). 
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Figure 16: The time dependence of the absolute value of e(t): y = \s(t)\ in 
the interval x G (1, 103). 

6 Verification of the relation HL{t)+Hs{t) — hu(t)+ 

All results in this Section have been rounded to the decimal third place. 
In accordance with formulae ( (1471) . (1451 ) for t = r L we have 

H L (t l ) ~ (497.648 - i ■ 4.458 x lO^MeV, 

fi s (r L ) ~ (497.648 - z • 2.471 x l(T 13 )MeV 

and the corresponding matrix elements of the effective Hamiltonian (formulae 
(I30|) -(j33l)) can be written for t = t l as 

h n (r L ) ~ (497.648 - i ■ 2 All x l(T 13 )MeV, (57) 
Mtl) ~ (1.787 x 10" 23 - £ • 6.401 x l()- 24 )MeV, (58) 
M T x) ~ (-1.799 x 10" 23 - i ■ 5.127 x 10~ 24 )iWW, (59) 
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h 22 (r L ) ~ (497.648 - i ■ 4.458 x 10~ 13 )Me\/. (60) 
For t = r L we get 

Mtl) + ^s(tl) = h u (r L ) + /i 22 (r L ) - (995.296 - % ■ 6.929 x l(r 13 )Mey.(61) 
Relation (1461) is also fulfilled at £ = r 5 

Hl(t s ) + Hs(r s ) = hn^s) + h 22 {r s ) ~ (995.296 - 2 ■ 9.623 x 10- 14 )MeV.(62) 



Comparing ( JB"Tj) with (l62l) we can see that 

Ml(tl) + Vs{tl) 7^ M r s) + /%(ts). (63) 
It is interesting to notice that 

Vl(t l ) + ii s (t s ) ^ (995.296 - i ■ 5.234 x 10~ 13 )MeV, (64) 

fj, L (r s ) + Hs(tl) ^ (995.296 - i ■ 2.658 x 10~ 13 )MeV, (65) 

hn(r s ) + h 22 {r L ) ~ (995.296 + i ■ 4.644 x 10~ 13 )MeV, (66) 

hu{r L ) + h 22 (r s ) ~ (995.296 + i ■ 3.247 x l(r 13 )MeV, (67) 

>From our calculations it follows that the real parts of formulae (1641) - (1671) 
differ in the twelfth or thirteenth decimal place. This means that 

Hl(ts) + I^s{tl) ^ ^l(t s ) + hs(t 3 ), (68) 

Pl[ts) + ^s{tl) ^ Vl(t l ) + fx s (r L ), (69) 

PlM + V>s{tl) ^ V>l{t l ) + HsM, (70) 

V>l{ts) + I^s{tl) ^ hn(r s ) + h 22 (r s ), (71) 

Vl(t l ) + ^s(rs) ^ h n (r L ) + h 22 (r L ), (72) 

A*l(tl) + IJ>s(ts) ^ hi(r L ) + h 22 (r s ) (73) 

and so on. 
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7 Final remarks 



First, as it was pointed out in [12], let us notice that in the considered 
model some relations assumed there and allowing to perform integrations 
of type ( fT7l ) are not valid in the K° — K° system (see a comment between 
formulae (5.10) and (5.11) of [12], Sec. 5). These relations were also used 
in [12] and |15| . Next, a drawback of our model is that at t = we obtain 
hn(t = 0) = oo and Ji22(t = 0) = oo. However, this model allows to study 
all the consequences of Khalfin's Theorem and theorems considered in [16j. 
Bearing in mind the limitations of our model mentioned above one should 
not expect that our calculations based on this model will result in an exact 
reconstruction of all experimental parameters characterizing the neutral K 
meson system. 

Results presented in Sec. 3 show how Khalfin's Theorem works. We can 
conclude that the effect of this Theorem should be visible in experiments 
with the neutral kaon complex in which the modulus of r(t), f TTTh . can be 
measured with the accuracy of order 10~ 16 or better. On the other hand, 
these results are in perfect agreement with the supposition formulated in 
[19] . Experimental results give |1 — \r(t)\ | ~ 1CT 3 = const with some limited 
accuracy (see, e.g. [20]) . The explanation of this fact proposed in [19] is 
based on the assumption that 

r{t) = r LOY + d(t), 

where r LOY stands for r{t) calculated within the LOY theory, and d{t) is 
assumed to be a function varying in time t such that \d{t)\ < 10~ n . 

Results obtained in Section 4 suggest that the real part of the difference 
(/in(£) — fi22{t)) is different from zero for very large times t: from t ~ 0, Itl 
up to t ~ 100t£. Moreover, after division by m average , this difference is 
only a little smaller than the corresponding experimental value [20]. The 
imaginary part of (hu(t) — h 2 2(t)) turned out to be different from zero as 
well. However, this part oscillates about 2 • 10 -13 MeV very fast. Note 
that from the results contained in [16] it follows that these differences should 
differ from zero for all £ > 0. Within the standard treatment of the neutral K 
system the measurement of the difference of masses (m K o— m^-o) is considered 
as the CPT invariance test. This interpretation of such tests is based on 
the properties (ED, © of the LOY approach: Within the LOY theory CPT 
symmetry is conserved only if (m^-o — m^o) = 0. The results obtained in 
Sec. 4 and in [15], [£6] [[7] show that such an interpretation of this test is true 
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only for the LOY approximation and beyond LOY approximation properties 
©, © do not occur. It seems to be obvious that the description of neutral 
K complex using the more accurate formalism than the LOY approximation 
leads to a more realistic description of such a complex. So, if within the 
more accurate theory one obtains (m^-o — ra^o) 7^ when CPT symmetry 
holds and CP symmetry is violated then one is forced to conclude that such a 
property must be valid in the real CPT invariant systems. Therefore, taking 
into account results obtained in Sec. 4 (and in [15], CEEl [17]) the conclusion 
that the measurement of the mass difference, [m K o — m^o), should not be 
considered as CPT invariance test seems to be correct. 

Results in Section 5 show that the imaginary parts of parameters fiL,^s 
and €l,£s vary in time too. We can say the same about their real parts. The 
oscillation amplitude is of the order of 1CT 13 for Im {^l(s)) and it is smaller 
than 10~ 27 for Im {eus\). The real parts Re (hl{s)) an d Re (el(s)) oscillate in 
a similar way as their imaginary parts and this is the reason why they cannot 
be shown in our figures. The parameter e, ( |54l) . is also a quantity varying in 
time as we can see from the graphical results (see FigH5] and FigJT6l). The 
absolute value of |e| oscillates around the value 8 • 10~ 6 . 

In Section 6, relation (1461) was investigated. We know, from general con- 
siderations, that it has to be fulfilled for every time t irrespective of whether 
we consider an approximate model of neutral K meson system or if we investi- 
gate the exactly solvable model of this system. The considered mathematical 
results show that the left side of equation f j46"l) is the same as its right side for 
t — Ti, t — Ts- A similar mathematical result was obtained for other times. 
From relation (l63l) it follows that the left side of equation (1461) , as well as its 
rights side, are not constant in the time. 

Relations (1641) - (1731 ) show that equation ( 1461 ) is no longer fulfilled when 
the separate components of sums appearing in its left and right sides are 
taken at different moments of time and then inserted into this equation. 
From this observation and from results in Section 4, we can draw an impor- 
tant conclusion concerning the methods of experimental data registering and 
experimental data processing. 

Let us note that an experimental system containing detectors which reg- 
ister the neutral K meson decay products can be schematically presented as 
in FigHTl This Figure presents a longitudinal section of a cylindrical vacuum 
chamber. A K meson stream is fired into this chamber along a horizontal 
axis / on the left side. Detectors D surrounding the chamber form its walls. 
These detectors register the neutral K meson decay products. In the first 
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region of this chamber (I), we observe a great amount of the decay products 
of neutral kaons into two pions (Kg — ► 2tt). In the second region of this 
chamber (II), we usually observe a great amount of the decay products of 
neutral kaons into three pions (Kl — ► 3tt). We can interpret the I axis 
as a path in an uniform straight-line motion of neutral kaons, whose decay 
products are registered by detectors D. We can write // = vk s -t~s in the first 
region (I) and In = vk l • t l in the second region (II) (where vk S{l) is the kaon 
Ks(L) speed). One can also obtain (for comparison): lj — c-Ts — 0.026805 m 
and l n — c ■ Tl — 15.51 m (c = 3 x 10 8 s - the speed of light in vacuum, 
r s = 0.8935 x 10~ 10 s a r L = 5.17 x 10~ 8 s). 
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Figure 17: A scheme of the experimental set for the experiment with the 
neutral K meson described in this section. 

Let us now return to the above mentioned conclusion. From (l64l) - ( l73l ) 
and from results presented in Section 5, it follows that only these parameters 
can correctly reflect real properties of neutral K system which are calcu- 
lated using only data obtained from a ring of detectors limited by distances 
(I, I + Al). Since t ~ I, the events are registered between t, t + At from 
the initial instant. Of course, Al should be as small as possible. In other 
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words, one should not use the experimental data obtained from the regis- 
tration of the neutral K meson decay products in the calculations if these 
neutral K meson decay products come from different and distant parts of the 
measurement set of the type shown in Fig{17l For example, one should not 
use calculations which were registered by the detector D in the region where 
t ~ t s simultaneously with the data which were registered by the detector 
D in the region where t ~ tl- One should not mix the data coming from 
different and distant parts of the measurement. If this rule is not observed, it 
may turn out that the obtained values fiL(s) (obtained on the basis of param- 
eters measured in this way) will not satisfy the consistency check given by 
( 1461) . Of course, in order to check the consistency of the experimental results 
with f T4"6l) the experiment should be conducted in such a way that both sides 
of (1461) can be found in independent measurements. Then one will obtain 
independently each other the left side and the right side of this equality. 
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